Antiferromagnetic spin-S chains with exactly dimerized ground states 
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We show that spin S Heisenberg spin chains with an additional three-body interaction of the 
form (Si_i ■ Si)(Si • Si+i) + h.c. possess fuUy dimerized ground states if the ratio of the three-body 
interaction to the biUnear one is equal to 1/(45(5' -I- 1) — 2). This result generaUzes the Majumdar- 
Ghosh point of the Ji — J2 chain, to which the present model reduces for S — 1/2. For 5 = 1, we 
use the density matrix renormalization group method (DMRG) to show that the transition between 
the Haldane and the dimerized phases is continuous with central charge c = 3/2. Finally, we show 
that such a three-body interaction appears naturally in a strong-coupling expansion of the Hubbard 
model, and we discuss the consequences for the dimerization of actual antiferromagnetic chains. 

PACS numbers: 75.10. Jm,75.10.Pq,75. 40. Mg 



Introduction ~ Over the years, exact results have 
proved to be extremely useful in quantum and statis- 
tical physics [1, 2]. In quantum magnetism, the Bethe 
ansatz solution of the spin- 1/2 Heisenberg chain [3] has 
led to the first proof that the spectrum is gapless [4], 
and its extensions, e.g., to the S* = 1 chain with bilinear 
and biquadratic interactions (BLBQ) with equal [5-7] or 
opposite [8, 9] amphtudes has helped a lot to clarify the 
physics of that model. In quantum frustrated magnetism 
[10], cases where an exact expression for the ground state 
wave function can be obtained have also played a very 
important role. For instance, for the spin-1 Heisenberg 
chain, the exact ground state of the AKLT point [11] has 
been a milestone in the confirmation of Haldane's predic- 
tion that the spectrum of integer-S' spin chains is gapped 
[12]. For spin-1/2 magnets, the first example of a gapped 
spectrum goes back to the Majumdar- Ghosh [13] (MG) 
point J2/J1 = 1/2 of the Ji — J2 model defined by the 
Hamiltonian 

nj,-j, = ^(Ji s, • s,+i + J2 s, • s,+2). (1) 

i 

At that point, the two fully dimerized states obtained as 
products of singlets on consecutive dimers and defined by 
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where \S{i, i + denotes the singlet formed by the spins 
at sites i and i + 1, have been shown by Majumdar and 
Ghosh to be exact ground states. Building on this result, 
it has been shown that the spectrum is gapped, and that 
this point is representative of an extended phase that cov- 
ers the parameter range 0.2411 < J2/J1 < +00 [14-16]. 
This seminal result has been at the origin of a long se- 
ries of experimental investigations of frustrated spin-1/2 
chains which started about 20 years ago with CuGeOa 
and which remains a very active field of research [17]. 



Attempts at generalizing the MG point to come up 
with a realistic model with fully dimerized states as ex- 
act ground states for larger spins have failed so far. The 
simplest idea is to consider the model of Eq. (1) for spins 
S >1 [18]. It is easy to convince oneself that the dimer- 
ized states of Eq. (2) remain exact eigenstates for any 
spin when J2/J1 = 1/2, but for 5* > 1, they are no longer 
the ground state. The problem can be traced back to the 
properties of a single triangle, into which the Hamilto- 
nian of Eq. (1) can be decomposed for J2/J1 = 1/2: For 
S = 1/2, the product of a singlet built out of two spins 
times any state of the third spin is a ground state. For 
5 > 1, the same state has a total spin 5, and it is not 
the ground state, which has total spin or 1/2 for integer 
and half-integer S respectively. 

Following Klein [19], an interesting alternative consists 
in building Hamiltonians as sums of local projectors on 
three spins to ensure that the product of a singlet with 
a single spin state is a local ground state. The simplest 
Hamiltonian of that kind takes the form [20] 
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where Pg^'^lf^'^ is the projector on the subspace of total 
spin S [33]. This projector can be written as 
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where the product runs from or 1/2 for integer or half- 
integer spins to 35". For S = 1/2, this Hamiltonian re- 
duces to the MG point of the Ji — J2 chain, but for 
S > 1, it is a polynomial in scalar products of pairs 
of spins of degree 35* or 3S — 1/2 for integer or half- 
integer spins, hence a very complicated Hamiltonian that 
seems difficult to realize in actual systems. The same re- 
mark applies to a spin-3/2 model recently investigated 
by Rachel [21], whose ground states are partially dimer- 
ized valence bond solid states, or to the generalizations 
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proposed by Rachel and Greiter[22] that lead to exactly 
trimerized resp. tetramerized ground states for S — 1 
resp. S ^ 3/2 models. 

In this Letter, we propose another generalization to 
arbitrary S of the spin-1/2 Ji — J2 model defined by the 
Hamiltonian 

H = ^ ( Ji S, • S,+i + J3 [(S,_i • S,) (S, • S,+i) + h.c.]) 

i 

(5) 

with Ji > 0. The number of sites N is assumed to be 
even, and we concentrate on periodic boundary condi- 
tions [34]. As we shall see, this Hamiltonian possesses 
for any value of S the equivalent of a MG point when 
J3/J1 = 1/(45(5 + 1) -2), at which the states ofEq. (2) 
are exact ground states, and it is realistic in the sense 
that it appears to next-to-leading order in the 1/U ex- 
pansion of the two-band Hubbard model that leads to 
the 5 = 1 Heisenberg model. 

For 5 = 1/2, it is easy to check that the Hamiltonian 
of Eq. (-5) reduces to that of Eq. (1) with J2 ~ J3/2. 
For 5 > 1, the three-spin interaction does not reduce 
to a next-nearest neighbor two-spin interaction, and the 
proof that the states of Eq. (2) are exact eigenstates is 
not a trivial extension of the MG proof. 

As in the 5 = 1/2 case, let us first determine un- 
der which condition the states of Eq. (2) might be exact 
eigenstates of Eq. (5). To be specific, let us consider 
\ipodd)- For i odd. Si • Si+ijV'odd) = -5(5 -I- l)\ipodd)- By 
contrast, for i even, the singlets on bonds {i — and 
{i + l,i + 2) are affected by Si • Si+i. However, the result- 
ing wave function does not contain states with arbitrary 
spin for the pairs {i — and {i + l,i + 2), but only 
triplets. Indeed, for two spins Si and S2, 5f|5(l,2)) 
is a triplet for all spin components a = x,y, z. This is 
clear for the z component since the SU(2) commutation 
relations imply that 

(5r + 52-)25r|5(l,2)) = 0, 
(5+ + 52+)(5r+52-)5f|5(l,2)) - 25J|5(1,2)), 

and by rotational symmetry, this has to be true of the 
other components as well. So, for i even, one can write 

Si ■ Si+il^orfd) = 

C„,ATa(i - 1, i)) \Ta'{^ + 1, * + 2)) [] '|5(j, J + 1)) 

cTjCT' j odd 

where the product over j is limited toj 7^ i — + and 
where the indices a, a' — 0,±1 keep track of the three 
possible triplets of a pair of spins. Since the total wave- 
function is a singlet, all coefficients must be equal to zero 
except Ci^-i, C_i_i and Cq.o, which must be related by 
Ci,_i = C_i_i = —C'ofi- Their common absolute value 
can be derived with the help of Clebsch-Gordan coeffi- 
cients, but this is unimportant for our present purpose. 
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Figure 1: (Color online) Spectrum of the Hamiltonian Hi 
[Eq. (7)] on three adjacent sites of the chain as a function of 
5(5'+ 1). The green line indicates the energy of the dimerized 
eigenstate on this three site system, E = —S{S + l)/2. 

The only relevant fact is that, since only triplets are in- 
volved, acting with Si_i-Si or Sj+rS,+2 on Si-S^+i [lAodd) 
will just multiply it by 1 — 5(5 + 1). This leads to: 

n\^odd) ^ ~^s{s + i)\tPodd) 

+ (Ji -(45(5+1)- 2) J3) J2 S,-S,+i|^odd)(6) 

i even 

If J3/J1 = 1/(45(5 + 1) - 2), the second term drops, 
and lipodd) is an eigenstate of H with energy per site 
— Ji5(5 + l)/2. Since the Hamiltonian is translationally 
invariant, this is also true for \ipeven)- 

To prove that these states are the ground states, let us 
decompose the Hamiltonian as H — Ji J2i with 

1 

+ 45(5+1) -2 ■ ■ + ^-"^-^ ■ 

The spectrum of this three-spin Hamiltonian can be 
worked out analytically for 5 = 1 and numerically for 
larger spin, with the result that the ground state energy 
EasCHi) is equal to -5(5+l)/2 (see Fig. 1). By the vari- 
ational principle, (n) > JlJ2^ EgsCHi) = ~NJiS{S + 
l)/2, a lower bound saturated by [V'odd) and \4'even)- This 
completes the proof that they are ground states of the 
Hamihonian of Eq. (5) when J3/J1 = 1/(45(5 + 1) - 2). 

Finally, it is plausible that these are the only ground 
states since the only ground states of Hi are the wave- 
functions with a singlet |5(i — 1,*)) or |5(i,«+l)), and the 
only common eigenstates are given by \ipodd) and \ipeven)- 
However, a mathematically rigorous proof that these are 
the only ground states for infinite systems would require 
an analysis similar to that of Ref.[ll] for the MG point 
of the spin-1/2 Ji — J2 model. 

Vicinity of the MG point for 5=1- We now con- 
centrate on the 5=1 model. At the Heisenberg point 
J3 = 0, the system is in the Haldane phase, which is 
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Figure 2: (Color online) Dimerization as a function of Js/Ji 
for different system sizes up to L = 250 sites in the vicinity of 
the phase transition J3/J1 ~ 0.11. The insets show the size 
dependence at J3/J1 — 0.1, 0.111 and 0.12, respectively. 



Figure 3: (Color online) Correlation length as a function of 
J3/J1 for different system sizes. Inset: spin and quadrupolar 
correlation functions of Eqs. (9) at the critical point J3/ Ji = 
0.111. 



gapped but not dimerized. Therefore, a phase transition 
has to appear between the MG point J3/J1 = 1/6 and 
the Heisenberg point. Let us investigate the nature of 
this transition numerically using the DMRG [23, 24]. 

The natural order parameter of this transition is the 
dimerization operator defined hy d — |(Si • S^+i — Si • 
Sj_i)| where + 1) is the central bond. Results for 
sizes up to 250 sites[35] are shown in Fig. 2. At the MG 
point, d is exactly equal to 2 for all sizes. The dimeriza- 
tion develops around J3/J1 = 0.11 in a way typical of a 
continuous transition. Assuming this to be the case, we 
have performed a finite-size scaling in the vicinity of the 
critical point, and we have identified the point where the 
dimerization decays to zero algebraically. This occurs at 
J3/J1 = 0.111(1) (middle panel of the inset of Fig. 2). 

This is further corroborated by our results for the cor- 
relation length, which we have obtained by fitting the ex- 
ponential decay of the spin-spin correlation function with 
x"^/^ exp (— x/^). The results up to 250 sites shown in 
Fig. 3 are consistent with a divergence at J3/J1 « 0.11. 
At the MG point, the correlation vanishes rigorously for 
all sizes. Together with the results for the dimerization, 
we therefore conclude that the transition is located at 
J3/J1 w 0.111. In the Supplemental Material, we also 
report on a scaling analysis of the fidelity susceptibility 
[1] that agrees with this estimate. 

Let us now try to further characterize the universality 
class of this phase transition. To this end, we have com- 
puted the central charge c from the block entropy of the 
system. Si = — Trgi^lng)^, with qi the reduced density 
matrix of a subsystem of size I. For a gapless system. 



Si = - In 
3 



'L . 




— sni 




TT 





5PBC, 



(8) 



in the presence of periodic boundary conditions, so that 
the central charge c is obtained by fitting the numeri- 



cal results to Eq. (11) [9]. Results for 50, 80 and 100 
sites [36] are shown in Fig. 4. They point rather convinc- 
ingly to c = 3/2. This suggests that the transition might 
be in the SU(2)fc=2 WZWN universality class [27], as 
the Takhtajan-Babujian (TB) point of the S' = 1 BLBQ 
chain, at which a transition from a gapped Haldane phase 
to a gapped dimerized phase takes place[8, 9]. 

To further test this conclusion, we have attempted to 
determine the scaling dimensions at the critical point 
which determine the exponents of the algebraic decay of 
the spin and quadrupolar correlation functions [28, 29] 

Cs{t,j) ^ {s^s;} ^ i^iy-^ , (9) 

Cq(z, j) ^ {S+f {Sjf + h.c.) ^ {^-])-^-'"\ 

For the SU(2)a;=2 WZWN transition, a = V^, i.e., the 
correlation functions decay with exponents 3/4 and 2, 
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Figure 4: (Color online) Fit of the Calabrese and Cardy for- 
mula [Eq. (11), continuous line] to the DMRG results (dots) 
for the block entropy at J3/J1 = 0.111. 
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Figure 5: Phase diagram of the Ji — J3 chain of Eq. (5) for 
5 = 1. 

respectively. A fit to the DMRG data at J3/J1 =0.111 
leads to exponents 0.72 and 1.83 for the corresponding 
correlation functions (see inset of Fig. 3), in reasonable 
agreement with the field theory prediction. Furthermore, 
the finite-size scaling of the correlation length at the crit- 
ical point is linear to a very good accuracy, which indi- 
cates that u = \. Finally, at the critical point we find 
d oc L~°-^^, implying fi/v ~ 0.47, hence /3 ~ 0.47 since 
V = \. In a related model, Nersesyan and Tsvelik [30] 
have predicted that the dimerization order parameter can 
be described as the product of four Ising fields. Three of 
them are ordered in the dimerized phase, one is disor- 
dered, and they are all critical at the transition point. 
Since the Ising exponent (3 is equal to 1/8, we expect the 
product of four critical Ising fields to scale with exponent 
j3 = 1/2. Again, the numerical estimate is in reasonable 
agreement with this prediction [37]. 

We therefore safely conclude that the MG point is rep- 
resentative of an extended phase which is separated from 
the Haldane phase by a continuous phase transition at 
J3/J1 ~ 0.111, and which extends to large values of J3, 
as in the 5 = 1/2 case [16]. The results are summarized 
in the phase diagram of Fig. 5. 

Discussion - Finally, let us discuss the implications 
of the present results for actual spin chains. For simplic- 
ity, we concentrate on spin-1 chains. [38] Starting from a 
two-orbital Hubbard model with repulsion U and Hund's 
rule coupling, a strong coupling expansion leads, to sec- 
ond order in the hopping integrals, to the S = 1 Heisen- 
berg model with bilinear coupling Ji . At fourth order, 
three extra terms appear: the three-body interaction J3 
of Eq. (5), a next-nearest neighbor bilinear coupling J2, 
as in Eq. (1), and a biquadratic interaction Juqi^i-^i+iY 
(see Supplementary Material) . The nature of the phase 
induced by these terms will depend on the microscopic 
parameters, but a reasonable case in favor of a spon- 
taneous dimerization in a realistic parameter range can 
be articulated around four points: 1) the J3 coupling 
generated by the fourth order perturbation theory is es- 
sentially always positive; 2) the critical ratio for dimer- 
ization J3/J1 = 0.111 is quite small and can be reached 
for reasonable values of t/; 3) the biquadratic interaction 
may be positive or negative. If it is negative, it favors 
dimerization. If it is positive, it is typically of the same 
order as J3, and preliminary results show that it should 
be significantly larger than J3 to suppress dimerization; 
4) to fourth-order, the next-nearest neighbor interaction 



is essentially ferromagnetic, and this would compete with 
dimerization. However, in actual antiferromagnets, it is 
in fact more likely to be antiferromagnetic due to residual 
direct superexchange, hence to be compatible with dimer- 
ization. So, we believe that the dimerization mechanism 
described by the model of Eq. (5) is a realistic potential 
source of dimerization in actual antiferromagnetic spin 
chains. We also note that for systems of ultracold alka- 
line earth atoms on optical lattices, higher order pertur- 
bation theory leads to the three-body term of Eq. (5) as 
well [31, 32]. In actual systems, this dimerization should 
be observable provided the interchain coupling and the 
temperature are both smaller than an energy scale of the 
order of the gap, a reasonable condition since the gap 
at the Majumdar- Ghosh point is expected to be a sig- 
nificant fraction of Ji (see Supplementary Material for a 
detailed discussion). 

Conclusions - We have shown that it is possible to 
generalize the spin-1/2 Ji — J2 model to larger spins in 
such a way that a Majumdar-Ghosh point where dimer- 
ized states are exact ground states is still present without 
making the model unrealistically complicated. For spin 
1, the additional interaction is a three-site term that ap- 
pears naturally at fourth order in a 1/U expansion of a 
two-band Hubbard model, and we have also shown that 
the MG point is representative of an extended dimerized 
phase separated from the Haldane phase by a continuous 
transition in the SU(2)fe=2 WZWN universality class. We 
hope that this new model will motivate the search for 
experimental realizations in quantum magnets and cold 
atoms. 
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The results can be extended to open boundary condi- 
tions, but then there is only one exact dimerized ground 
state. 

The Dimerization was computed with OBC. Close to the 
transition, we keep up to 1400 states and performed 9 
sweeps. This keeps the discarded weight in the last sweep 
below 10"" 

The central charge was computed with PBC. We keep up 
to 2500 states and performed 24 sweeps. This keeps the 
discarded weight in the last sweep below lO"'^ 
Note that the same procedure applied at the TB point 
of the S — 1 BLBQ chain yields a similar value for the 
exponent (0.46). 

Similar arguments should apply to arbitrary S since the 
fourth-order terms of the strong coupling expansion have 
to be the same by symmetry and counting arguments. 



Supplemental material to "Antiferromagnetic 
spin-5 chains with exactly dimerized ground 
states" 



FIDELITY SUSCEPTIBILITY FOR THE S = 1 
Ji - Js CHAIN IN THE VICINITY OF THE 
PHASE TRANSITION 



In this section, we discuss in more detail our findings 
for the fidelity susceptibility [1] ( Ji = 1 in the following) 



1-\{MJ3)\MJ3 + SJ3))\ 



(10) 



with \ipo{Jz)) the ground state of the Ji — J3 chain. Ac- 
cording to the analysis of Refs. 2, 3 and the numerical 
findings of Refs. 4-8, in the thermodynamic limit x{Jz) 
should either possess a divergence or a minimum at the 
critical point, depending on the values of the scaling di- 
mensions and of the critical exponents. The results of 
Fig. 6 indicate that in the present case a peak develops. 
We perform an extrapolation of the peak with system 
size and find that in the thermodynamic limit indeed a 
divergence is obtained. As can be seen in the inset of 
Fig. 6, extrapolating the peak position leads to a value 
Jc ~ 0.11, in agreement with the findings for the dimer- 
ization and for the correlation length. 



>^20 




Figure 6: (Color online) Fidelity susceptibility [Eq. (10)] as 
a function of J3/J1 for different system sizes in the vicinity 
of the phase transition J3/J1 ~ 0.11. The inset shows the 
finite-size scaling of the peak position. 
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Figure 7: (Color online) Central charge obtained for systems 
of size L — 50, 80 and 100 from a fit to Eq. (11) after different 
numbers of sweeps of the DMRG algorithm. The error bars 
only take into account the precision of the fit for a given 
number of sweeps. As one can see from the evolution of the 
value of the central charge with the number of sweeps, the 
values for L = 80 and L = 100 have not fully converged 
due the truncation of the Hilbert space, and the values for 
the largest number of sweeps we could achieve are probably 
underestimates. Because of these uncertainties, a meaningful 
finite-size scaling is not possible. Still, these results are clearly 
compatible with a central charge c — 3/2. 



0.7 
0.6 
0.5 
0.4 

:o.3 

0.2 
0.1 
0.0, 









0.5 



1 

J:i/Ji 



1.5 



Figure 8; (Color online) Gap of the system as a function of 
Ji/Ji for different system sizes. 



up to 144 sites. Moreover, one should keep in mind that 
the critical ratio J^/Ji is not known exactly, another po- 
tential source of error. Still, the results clearly point to 
a central charge c = 3/2. 



FINITE SIZE RESULTS FOR THE CENTRAL 
CHARGE 

To characterize the transition, we have computed the 
central charge c from the block entropy of the system. 
Si — — Trpf Inpf, with the reduced density matrix of 
a subsystem of size t. For a gapless ID system with 
periodic boundary conditions, it behaves as [9] 



Si 



In 







— sm 
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(11) 



Due to the limitations of the DMRG when dealing with 
systems with PBC, we are restricted to system sizes of 
the order of 100 lattice sites. In Fig. 7, we show the re- 
sults for systems with L — 50, 80 and 100 lattice sites 
including the error bars that come from the fit to the 
Calabrese-Cardy formula of Eq. (11). Another source 
of error comes from the truncation of the Hilbert space 
in the DMRG algorithm. To give the reader an idea of 
this error, we have plotted the central charge obtained 
after different numbers of DMRG sweeps up to the max- 
imal number of sweeps we could achieve. On the very 
fine scale of the plot, the change in the central charge 
is negligible for 50 sites, but it is already noticeable for 
80 sites and quite significant for 100 sites. Given these 
uncertainties, and the smallness of the deviations from 
c=3/2, it does not appear meaningful to perform a finite 
size extrapolation. This sould be contrasted to the spin- 
1/2 case treated in [10], where a finite-size analysis could 
be performed thanks to the good convergence achieved 



GAP AND CONDITION ON INTERCHAIN 
COUPLING AND ON TEMPERATURE 



At the Majumdar- Ghosh point, as in the case of the 
Ji — J2 model, the gap is expected to be a significant 
fraction of Ji [11]. We have checked this expectation 
with DMRG in the S = 1 case. Indeed, as shown in Fig. 
8, the gap increases very fast above J3/ Ji = 0.11 to reach 
values of the order of 0.7Ji around the MG point, a bove 
which it decreases to stabilize around 0.3-\/ + J^. 

In antifcrromagnets, the conditions to observe the 
dimerization are thus expected to be met in sufficiently 
anisotropic systems and provided one can reach tempera- 
tures smaller than the dominant coupling constant, which 
is essentially always the case. 

In cold atoms, to suppress the interchain coupling, one 
should work in a ID trap. The condition on the temper- 
ature should be translated into a condition on entropy. 
It is not possible to be quantitative without calculating 
the temperature dependence of the entropy, a calculation 
far beyond the scope of the present analysis, but to reach 
a temperature equal to a fraction of the main coupling 
constant means to work with an entropy per site equal to 
a fraction of ln(25' -I- 1). This is the typical condition to 
observe antiferromagnetic correlations in cold-atom real- 
izations of the Heisenberg model, and this is an issue on 
which the experimental community is currently actively 
working. 
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DERIVATION OF THE Ji - Ja MODEL FROM A 
TWO-ORBITAL HUBBARD MODEL 

The present section is devoted to the derivation of 
the effective Ji — J3 model from a microscopic Hubbard 
model. In solid state physics spin-1 systems may be real- 
ized in the case of transition metal compounds: If there 
are 2 degenerate orbitals due to the crystal field splitting 
and if the Hund's coupling Jh and the Coulomb repul- 
sion are large enough, at half-filling (i.e. two spin-1/2 
per transition metal ion) the system consists of localized 
spin S=l moments. 

Assuming this to be the case, we derive an effective 
spin S=l Hamiltonian up to fourth order in degenerate 
perturbation theory in the strong coupling limit of a two- 
orbital Hubbard model on a chain, and we show under 
which conditions on the original microscopic parameters 
(hopping integrals. Coulomb repulsion, Hund's coupling) 
this effective Hamiltonian reduces to a Ji — J3 model 
with nearest-neighbor coupling and three-body interac- 
tion terms. 



where S are spin-1 operators. The sum over < i,j > 
runs over nearest-neighbor pairs, the one over i,j ^ 
runs over next-nearest neighbor pairs, while the one over 
< i, J, A: > runs over all sequences of three spins. 



The effective Hamiltonian consists of different terms. 

At lowest order, one recovers the Heisenberg model with 

(2) 

a nearest-neighbor spin coupling Jfj(,js which gets renor- 

malized at fourth order by the coefficients ^^oL- 
fourth order, two other 2-site terms appear: A next- 
nearest neighbor coupling Jnn\ as in the S=^ case [14], 

and a biquadratic coupling Jg^qu typical oi S — 1 sys- 
tems. Finally, there is in additional three-site interaction 
which cannot be reformulated as a 2-site operator. 
Let us mention that these terms, which appear in the 
perturbation expansion of the Hubbard model, have also 
been extracted from ab initio calculations in a different 
context [15]. 



Generalized Hubbard model 



Our starting point is the following Hamiltonian at half- 
filling: 

'^"^ = H H ^t'i,m'4maC3m'a (12) 
i,j m,rn' o 

~t- 2 ^ ^ ^ ^ ^ ^ ^ram''^itna'^im' a' 
i ra,m' cr,CT' 

+ 2Jpcl^,^,cl^,^Cimi7'Cim<j} 

where i,j are the site indices, m,m' refer to the or- 
bitals a and 6, and a to the electronic spin. The hopping 
integrals between two neighboring orbitals are denoted 
by ^, , the on-site Coulomb repulsions by Umm' , Jh 
represents the Hund's coupling and Jp the pair hopping 
amplitude. Furthermore, we assume that additional re- 
lations, typical of cubic symmetry, are satisfied, namely 
Uaa — and U — Uaa ~ '^Jh- This kind of Hamiltonian 
has been extensively discussed in the context of systems 
with orbital degeneracy [12, 13]. 
Effective spin-model to fourth order in degenerate 
perturbation theory 



In order to have more compact expressions, we intro- 
duce the following relations: 



''aa ^ '' 



bb 



''aa''ab ^ ''bb'-ab 



,2 ,2 

''aa'-bb 



iaa^bb^ab 



Using degenerate perturbation theory, the effective ^he various coefficients of (13) can then be expressed 
spin model of Eq. (13) on a chain takes the form terms of the microscopic parameters of the original 

ff = ^(2) _|_ ^(4) Hubbard model as sums of terms classified according to 

jj{2) _ j(2) /g g N the combination of U, Jh and Jp that appears in the 

Heis Z-^ ^ * J' rlpnnminntnr' 



<t,J> 

H^'' = E (^^fis.-S,+J^JjS..S,)^ 

<i,j> 

+ E/ "^nni^i ' ^JJ 



7(4) 
'Heis 



23h + U 



32t| 



8t| 



{2Jh + U)^ 



+ 



2tf 



9^4 



+ 



16t| 



+ 



;(4) 
'Biqu 



+ 



7(4) 
■^3 



+ 



+ 



-8t|p + 16t| - 8tl 



ab 



+ 



-4tt 



- 8t| - At* 



ab 



+ 



-16t|-4t^p-8t|-4t^, 



(-5Jh - 2Jp 



2tt 



-8t| 



(7) (2 J;, 

4t4 



+ 



+ 8t! - 4t 



afc 



?7)2 (-5J,, + 2Jp-i7)(2Jft + C/)2 

-16t| - Wtj 



16t^ - 16t| 



(2J,. + [/)= 



(-3J,, - 2Jp)(2Jft + C/)2 (-5Ja - 2Jp)(2J^ + [/)- 



+ 



3t^ 
2 



2t|- 



5t: 



2p 

2 



''ab 



3h{23h + U)^ 

K> + 



+ 



8t4 



44 



(-6J^ + 4J,)(2J^ + C/)^ 



+ 



4t,;. 



4t 



2tf 



(-6J,,-4Jp)(2J;, 

4tt + 16t| + 24j, 

{23h + Uf 
t* 



u{23h + uy 



4t| 



+ Kb 



+ 



tt-t 



J/i(2J/i 



4t^p 



8t|- 



■4t4 



+ 



2tf 



■4t|- 



2p 

-C/)2 



8t| 



(-3J/, 
+ 



2Jp)(2J;, + ;7)- 



8t4. 



■4d 



+ 



2t^, 



{-23h-2,U){23h + UY {-A3h-U){23h + Uy (-5Jh - 2Jp - C/)(2J^ + C/)^ 
(-5Jfc + 2Jj,-[/)(2J^ + ?7)2 



+ 



-4tt - 16t| - 8t| - 4t^fc -2tf + 2t|p 16t| + 16t| 

(2j^ + c/)3 ^ jft(2j;, + uy ^ (-3J,, - 2jp)(2j,, + uy ^ 

-4ttp + 8tt-4t4, _2t4 + 4t|-2t4, -8ti - 2t|p - 4tt - 24 

(-2Jft - 3C/)(2J^ + Uy ^ (-4J^ - C/)(2J,, + Uy ^ (-5Jfc - 2Jp - f/)(2Jft + ?7)2 

-2t|^ + 4t|-24 
(-5J;, + 2Jp-[/)(2J^ + i7)2- 
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